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I. INTRODUCTION 



q3 , 

r-{ , Analyzing pure quantum states in an infinite dimensional Hilbert space it is useful to distinguish a family of coherent 

>—i states, localized in the classical phase space and minimizing the uncertainty principle. These quantum analogues of 

fvi points in the classical phase space are often considered as 'classical' states. For an arbitrary quantum state one may 

pose a natural question, to what extent it is 'nonclassical' in a sense that its properties differ from that of coherent 

T-H ' states. In other words, is there any parameter that may legitimately reflect the degree of nonclassicality of a given 

^ , quantum state? This question was motivated with the first observation of nonclassical features of electromagnetic 

'^ ' fields such as sub-poissonian statistics, antibunching and squeezing. Additionally, it is well known that the interaction 

of (non)linear devices with quantum states may flip from one state to another; for instance, nonlinear devices may 

produce nonclassical states from their interaction with the vacuum or a classical field. A systematic survey of 

^^ nonclassical properties of quantum states would be worthwhile because of the nowadays ever increasing number of 

^+ I experiments in nonlinear optics. An earlier attempt to sheding some light on the nonclassicality of a quantum state 

f^ ' was pioneered by Mandel [l|, who investigated the radiation fields and introduced a parameter q measuring the 

^H , deviation of the photon number statistics from the Poissonian distribution, characteristic of coherent states. 

Oj' In general, to define a measure of nonclassicality of quantum states one can follow several different approaches 

' ! 0- Distinguishing a certain set C of states (e.g. the set of coherent states |a)), one looks for the distance of an 

r"| ' analyzed pure state {ip) to this set, by minimizinga distance d(|V'), \ct)) over the entire set C. Such a scheme based 

Cd ' on the trace distance was first used by Hillery 0, |j| , while other distances (Hilbert-Schmidt distance _5^, _6j or Bures 

^ . distance |3;la|) were later used for this purpose. The same approach is also applicable to characterize mixed quantum 

^ ' states: minimizing the distance of the density p to the set of coherent states is related |a, to the search for the 

.5^ ' maximal fidelity (the Hilbert-Schmidt fidelity Tr (pa) or the Bures-Uhlmann fidelity ( Try^p^/^crp^/^ j ) with respect 

^ ' to any coherent state, a = |a)(Q;|. On the same footing, the Monge distance introduced in [i3,lll| may be apphed to 
describe, to what extent a given mixed state is close to the manifold of coherent states. 

Yet another way of proceeding is based on the generalized (Cahill) phase space representation Rr of a pure state, 
which interpolates between the Husimi (Q), the Wigner (W) and the Glauber-Sudarshan (P) representations. The 
Cahill parameter r is proportional to the variance of a Gaussian function one needs to convolute with P representation 
to obtain Rr [l3|. In particular for t = 1,1/2,0 one obtains the Q-, W- and P- representations, respectively. By 
construction the Q representation is non-negative for all states, while the Wigner function may admit also negative 
values, and the P representation may be singular or may not exist. 

The smoothing effect of Rr is enhanced as r increases. If r is large enough so that Rr becomes positive definite 
regular function, thus acceptable as a classical distribution functioii, then the smoothing is said to be complete. The 
greatest lower bound r^ for the critical value was adopted by Lee [Ij, [3 , as nonclassical depth of a quantum state 
and this approach was further developed in |l5l Il6l Il7| . The limiting value, Tm = 1, corresponds to the Q function 
which is always acceptable as a classical distribution function. The lowest value, r^ = 0, is ascribed to an arbitrary 
coherent state because its P function is a Dirac delta function, so its e-smoothing becomes regular. The range of r^ 
is thus T,„ e [0, 1]. 

If the Husimi function of a pure state admits at least one zero (5(ao) = 0, then a Cahill Rr distribution with a 
narrower smearing, r < 1, becomes negative in the vicinity of ao. Therefore the classical depth for such quantum 
states is maximal, r™ = 1 .1^]. The only class of states, for which Q representation has no zeros, are the squeezed 
vacuum states, for which r^ is a function of the squeezing parameter s. In the limiting case s — one obtains the 
coherent state for which the Rq = P distribution is a Dirac delta function, that is r™ = 0. 



A closely related approach to characterizing quantum states is based on properties of their Wigner functions in 
phase space {p,q}- One can prove that the Wigner function is bounded from below and from above [l^l- In the 
normalization // W{q,p)dqdp — 1 used later in this work, such a bound reads |VF((7,p)| < l/nh. Further bound on 
integrals of the Wigner function were derived in [ig , while an entropy approach to the Wigner function was developed 
in [l9. 20], 

In order to interpret the Wigner function as a classical probability distribution one needs to require that W is non- 
negative. As found by Hudson in 1974 [21|, this is the case for coherent or squeezed vacuum states only. A possible 
measure of nonclassicality may thus be based on the negativity of the Wigner function which may be interpreted as 
a signature of quantum interference. 

The negativity of the Wigner function has been linked to nonlocality, according to the Bell inequality (22j| , while 
investigating the original Einstein-Podolsky- Rosen (EPR) state |23- In fact Bell argued that the EPR state will 
not exhibit nonlocal effects because its Wigner function is everywhere positive, and as such will allow for a hidden 
variable description of correlations. However, it is now demonstrated |2J, |23 that the Wigner function of the EPR 
state, though positive definite, provides a direct evidence of nonlocality. This violation of the Bell's inequality holds 
true for the regularized EPR state 26| and also for a correlated two- mode quantum state of light 'T^ ■ 

It is also worth recalling that the Wigner function can be measured experimentally 28] , including the measurements 
of its negative values |29j. The interest put on such experiments has triggered a search for operational definitions of 
the Wigner functions, based on experimental setup |30ll3ll|. 

The aim of this letter is to study a simple indicator of the nonclassicality, which depends on the volume of the 
negative part of the Wigner function. To demonstrate a potential use of such an approach we investigate certain 
families of quantum states. The nonclassicality indicator is defined in section 2. The Schrodinger cat state, being 
constructed as coherent superposition of two Gaussian wave packets, is analyzed in section 3 while section 4 is 
devoted to Fock states and to the squeezed displaced Fock states. Finally in section 5, a brief discussion of results 
and perspectives is given. 

II. THE NONCLASSICALITY INDICATOR 

The Wigner function of a state \ip) defined by |33.l33l| 

1 f^°° X X 

Wi,{q,p)^— dx{q~-\ij){'^\q+-)ex.p{ipx) (2.1) 

satisfies the normalization condition JJ W^{q,p)dqdp = 1. Hence the doubled volume of the integrated negative part 
of the Wigner function may be written as 

<5(V) = J J [miq,p)\ - W^{q,p)] dqdp = JJ \W^{q,p)\dqdp - 1 . (2.2) 

By definition, the quantity S is equal to zero for coherent and squeezed vacuum states, for which W is non-negative. 
Hence in this work we shall treat S as a. parameter characterizing the properties of the state under consideration. 
Similar quantities related to the volume of the negative part of the Wigner function were used in [SJ, |33, |33 to 
describe the interference effects which determine the departure from classical behaviour. 

Furthermore, a closely related approach was recently advocated by Benedict and collaborators [33, ISg- Their 
measure of the nonclassicality of a state \ip) reads 

where /+ {ip) and /_ {ip) are the moduli of the integrals over those domains of the phase space where the Wigner function 
is positive and negative, respectively. The normalization condition implies /+ — /_ = 1, so that ly ~ 2/_/(2/_ + 1) 
leads to < :^ < 1. Using this notation we may rewrite l|2.2|l as J = /++/_— 1 = 2/_ and obtain a simple relation 
between both quantities 



2/_ 



1 + 2/_ 1 



(2.4) 



with 6 = 1^/(1 — ly). It turns out that both quantities are equivalent in the sense that they induce the same order in 
the space of pure states: the relation (S(V'i) > Siil'2) implies ^{ipi) > J^('02)- However, from a pragmatic point of view 
there exists an important difference between both quantities. 



To compute explicitly the quantity H2.3|l one faces a difRcult task to identify appropriately the domains, in which 
the integration has to be carried out. On the other hand, knowing the Wigner function W{q,p) of a quantum state, 
it is straightforward to get its absolute value and to evaluate numerically the integration (|2.2(l . 

Let us emphasize again that the Hilbert space containing all pure states is huge, so one should not expect to 
characterize the nonclassical features of a quantum state just by a single scalar quantity. Our approach focuses on a 
particular issue, whether the Wigner function is positive and may be interpreted as a classical probability distribution. 
Therefore, the proposed indicator S should be considered as a tool complementary to these worked out earlier and 
reviewed above. 

III. THE SCHRODINGER CAT 

A quantum state, called Schrodinger cat, is a coherent superposition of two coherent states localized in two distant 
points of the configuration space, ±go. The wave function of such a state reads in the position representation 

M^(g) = ^ [</-+(<?) + 0_(g)] (3.1) 



where 

0±(9) = (^) exp(- — (g±go)'+*Y(g±go)) ■ (3.2) 
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From now on atomic units are used (m = h = to = 1). In other words we measure the size of the product pq in units of 
h. The classical limit /i — > means the action pq characteristic of the system is many order of magnitude larger than 
h. A glance on Eq. (|3.2II reveals that the phase, governed by po, is of great importance in that it induces oscillations 
on the wave function as can be seen in Fig.l. Note that the normalization constant TV depends on the location of the 
centers {qo,po) of both coherent states that make up the cat state. Therefore one sees that the Wigner function may 
depend not only on the distance 2qo between the both states, but also on their momentum, pq. So far, the studies on 
the cat states J3l| were usually restricted to the case of standing cats, po = 0. In this letter we demonstrate that the 
parameter po influences the shape of the Wigner function, in particular, if go ~ 1 and both packets are not spatially 
separated. 

Inserting (|3.1|) into the Wigner function (|2.1|) one obtains 

W^iq,p) ^ W+iq,p) + W^{q,p) + W,nt{q,p) ■ (3.3) 

Here 

W±{q,p) ^—exp{~{q±qof-{p^Pof) (3.4) 

represent two peaks of the distribution centered at the classical phase space points (itgoiPo); while 

Wint (<?, p) = — COS (2mo) exp {-q^ - [p - Pof) (3.5) 

■K 

stands for the interference structure which appear between both peaks. Normalizing (|3.1|) yields 

iV-(l + cos(2pogo)exp(-g2))-i/2 _ (gg^ 

Making use of the formula H3.3|) for the Wigner function of the cat state |^) its nonclassicality parameter 

(5(*) = jj \W+(q,p) + W,,u{q.p) + W^(q,p)\dqdp - 1 (3.7) 

may be approximated by 



6{^) « iV' 



/7 
-j= |cos (2p(jo)| exp (-(p - pof 



1 . (3.8) 



Strictly speaking the right hand side of equation (|3.8|) forms an upper bound for i5(^), which may be practically 
used as its fair approximation. Because of the oscillations of the absolute value of cosine, it is difRcult to perform 




FIG. 1: Schrodinger cat states wave functions plotted with po — (left) and with po = 4 (right). Dashed and solid lines 
represent the imaginary and the real part of the wave function, respectively. Notice that the envelopes of both wave functions 
do coincide. 




FIG. 2: Plots of the Wigner functions of the Schrodinger cat states l|8..'i|l . Each panel is labeled by the separation distance qo, 
the momentum po and the resulting indicator 6. Observe that for intermediate separations, go ~ 1, the indicator 5 changes 
with Po- Upper row shows the 'standing cats' (po — 0) while the cats in motion (po — 4) are represented in the lower row 



the integration analytically. In the special case qo = 0, the superposition of coherent states H3.1() reduces to a single 
coherent state and correspondingly (|3.8|l leads to (5(^) = 0. 

Fig. 2 shows plots of the Wigner function of the cat states for several values of the separation qq and the momentum 
Pq. One clearly sees the formation of the quantum interference structure halfway between the two humps as the 
separation distance qo increases. The frequency of the interference structure increases with the separation 131]. For 
intermediate separations {0 < qo < 4), the Wigner function changes its structure with po, see fig. 2b and 2c. However, 
for a larger separation distance, qo > 4, the Wigner function for po = Pi ^ may be approximated by the Wigner 
function for the state with po = translated by a constant vector Ap — pi . 

In the case of 'standing cats', {po = 0), the indicator S increases monotonically with the separation qo, and reflects 
presence of the interference patterns at q = - see Fig. 4k. The growth of the nonclassicality saturates at go ~ 4, 
as the interference patterns become practically separated from both peaks, and the parameter S tends to the limiting 
value, Smax ~ 0.636. In the limit qo ^ oo the oscillations of the cosine term in Eq. 1)3. 8|l become rapid and a crude 
approximation | cos(goPo)| ~ 1 gives an explicit upper bound S < 2iV^ — 1 « 1. 



H^ 




FIG. 3: As in Fig. 2. The values of go labeling each panel, correspond to the successive extrema (g,h,i) of the indicator 5 
plotted in Fig. 4m as a function of go. Wigner funtion for the 'Cats in motion' as in lower row of Fig. 2, for the selected values 
of po, for which the dependence S{qo) achieves its extrema. 




FIG. 4: Indicator S of the Schrodinger cat state \xf}) as a function of the separation distance go and several values of po as 
labelled on each panel. Grey dots (a-f ) refer to labels of individual panels of Fig. 2. while grey dots (g,h,i) refer to that of fig. 3. 

This picture gets more complicated for the states with po ^ 0, in particular for a small separation distance, 
(0 < go ^ 4). In this case, 6 exhibits oscillations as shown in Fig. 41, 4m, 4n. To shed some light on this behavior 
we have chosen to plot in Fig. 3 the Wigner function for which 5{qo) achieves extremal values. For instance, S at 
qo = 0.725 (Fig. 3b) is smaller than at qo = 0.4 (Fig. 3a) or 1.175 (Fig. 3c). This is due to the the interference 
structure, which is not symmetric with respect to the reflection p — > —p, in contrast to the case of cats with pq = 0. 

As shown in Fig. 4, the frequency of oscillations increases with pq, but the limiting value S{qQ —^ cxi) does not 
depend on the initial momentum pQ. This can also be demonstrated, investigating the dependence of the quantity S 
as a function of po- As follows from Eq. (|3.8f) . the indicator S displays regular oscillations with the period Posc = tt/qo 
- See Fig. 5. In other words a non-zero separation parameter go breaks the translational invariance in momentum 
and introduces a characteristic momentum scale Posc ^ l/9o- Note that the amplitudes of the oscillations decrease 
fast with qo, so that for well separated cats with qo > A the quantity d is practically independent on pq. 



IV. GENERALIZED FOCK STATES 



Let us consider the squeezed displaced Fock state defined by 

\/3,r],n) ^ Sifi)D{p)\n) , 



(4.1) 



where \n) is the original Fock state and n = 0, 1, 2, . . . . The displacement D{f3) and the squeezed S{ri) operators are 
defined by 0,121 



D{(3) := exp(/3a^ - f3*a) and S{tj) := exp ( |(»7*a^ - W^) ) 



(4.2) 
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FIG. 5: Indicator 5 of the Schrodinger cat state j'J') as a function of the momentum po for certain values of the separation go- 

where a and a' are usual photon annihilation and creation operators, respectively. The complex variable P represents 
the magnitude and angle of the displacement. Similarly, writing the complex number in its polar form, i] = sexp(i0), 
it is easy to see that the radius s p lays the role of the squeezing strength while the angle (j) indicates the direction 
of squeezing. It was shown in [l2| that the displacement operators £>(/?) form a complete set of operators. Thus 
any bounded operator F, (for which the Hilbert-Schmidt norm || F ||— y''Tr{F'i F) is finite), can be expressed in 
the form F = J f{^)D^^{(^)(P^/TT in which the weight function /(^) = Tr(FZ?(^)) is unique and square-integrable. 
Given that every density operator is bounded {Ti{p' p) — Ti{p^) < 1), one may write an arbitrary density operator 
p = J x{C)D~^ {£)<P £, / TT ■ Here the weight function x(^) = Tr(/9Z?(^)) is just the expectation value of the displacement 
operator commonly known as characteristic function. The complex Fourier transform of x(f) defines the Wigner 
function 



W{a) 



/ — x(C)exp(ra-Ca* 



(4.3) 



One may therefore express xiO ^^ terms of the Wigner function by performing the inverse Fourier transform as 



xiO ^ J ^W{a)exp{^a* - Ca) , 



so that upon substitution into the density operator expression above, one gets 



-W{a)T{a) 



(4.4) 



(4.5) 



The operators T{a) denote 



= 2D{a){-~iy^''D-\a) , 



(4.6) 



so that the Wigner function may be interpreted as a weight function for the expansion of the density operator in terms 
of the operators T{a) L2j . These operators are Hermitian, T = T\ and possess the same completeness properties as 
the displacement operators D{a). Making use of the parity operator (—1)" " = exp(i7ra^^a), one finally shows that 



W{a) = 2(-l)"Tr(pi:>(2a)) 

with n = a^a being the photon number. 

In the case of the squeezed displaced Fock states, p ~ \P,ri,n){(3,ri,n\, the Wigner function becomes 



(4.7) 



W„ia) = 2(-l)"(/3,77,n|7^(2a)|/3,77,n) 



(4.8) 
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FIG. 6: The nonclassicality indicator 5{\n)) of the Fock states versus the quantum number n < 250 (sohd line). Dashed hue 
represents ^\/n plotted for comparison. 



Performing explicitly calculations of matrix elements, one obtains : 

W^{a) = -(-!)" exp(-2|6p)i„(4|6n 



(4.9) 



with b = cosh(s)(a* — /3*) + exp(— i0) sinh(s)(a — /3)[40|- Here L„ denotes the Laguerre polynomial of the 7i-th order. 
The Wigner function (|4.9I) allows us to compute the nonclassicality parameter 5{\(3,ri,n)) for a given displaced 
squeezed Fock state \/3,ri,n). In what follows certain special cases will be investigated such as squeezed displaced 
vacuum states, pure Fock states and squeezed displaced Fock states. It will be therefore convenient to represent the 
complex variable a by the position and momentum coordinates, a — ~7=('Z + w)^ ^-nd treat likewise the displacement 

operator, p = -^{qq + ipa). 



V2^ 

Substituting P — rj = in eq. 



Wniq,p) = 



yields the Wigner function for the Fock state \n), 

i-iy 



■ exp [- 



-p2)]i„[2(q2+p2)] 



(4.10) 



This allows to evaluate analytically the indicator S{\n)), for n = 1,2,3,4 

"^(lO)) — (vacuum) 

S{\1)) == -i- - 2 « 0.4261226 

e ' 

S{\2)) = 4((2 + \/2)e~'~^ + (-2 + y2)e"'+^) 

6i\3)) w 0.97667 
J(|4)) « 1.19138 , 



0.72899 



(4.11) 



since the zeros of the Laguerre polynomials are available up to the 4-th order. For larger n we computed the quantity 
6{\n)) numerically and plotted in Fig. 6. The indicator S grows monotonically with n, as the number of zeros of the 
Laguerre polynomial L„ increases with n. For n e [1,250] this dependence may be aproximated by ■^y/n. Hence, 
the larger the quantum number n, the less the Wigner function W|n) can be interpreted as a classical distribution 
function. 

Setting n = in (|4.9|) one obtains a squeezed coherent state or squeezed vacuum state. Choosing the squeezing 
angle (j) — 0, one sees that the Wigner function is a Gaussian centered at the displacement vector (qotPo) with the 
shape determined by the squeezing parameter s. 



Wo{q,p) = -exp 



^e^'{q~qo)' 



[P-Pof 



□ 2s 



(4.12) 



In such a case the Wigner function remains everywhere non-negative for any choice of the squeezing and displacement 
parameters |2lj, so that the nonclassicality indicator vanish, S{\P, s, 0)) = 0. Note that the displacement of any state 




FIG. 7: Contour plots of the Wigner functions of the squeezed Fock states |0, s exp(i7r/6), 3) labeled by the squeezing strengths 
s. Irrespectively of s the indicator 5 « 0.97. 

in phase space does not change the shape of the Wigner function, so the quantity 6 is independent of the displacement 
operator D{(3). 

Furthermore, the squeezing operator S{ri) influences the shape of the Wigner function, but does not lead to a change 
in the volume of its negative part. Therefore, the parameter 6 does not also depend on the squeezing. As an illustration 
we have chosen the squeezed (|a| — s,(j) — 7r/6) displaced (/? = 0) third photon (n = 3) state, |0, sexp(i7r/6), 3). The 
contour plots of the Wigner function of such a state are shown in Fig. 7 for some values of the squeezing parameter s. 
The indicator 6 is equal to 0.9762, irrespective to the squeezing strength. If squeezing is strong enough, the ring-like 
Wigner function collapses to a quasi one dimensional object with a cigar form. 

The squeezed vacuum is often described as a nonclassical state |3l| . Since the quantity 6 does not depend on 
squeezing, it should not be interpreted as the only parameter which characterizes the nonclassicality. To describe the 
nonclassical features of the squeezed states one may use, for instance, the nonclassical depth [isL llsL Il7| . 



CONCLUDING REMARKS 



In this work we have proposed a simple indicator of non-classicality which measures the volume of the negative 
part of the Wigner function. Although the proposed coefhcient (5 is a function of the related quantity ly, recently 
introduced by Benedict, Czirjak et al. 37. 38]. it is much easier to compute numerically. 

The quantity l|2.2|l was used to analyze exemplary quantum states, including the Schrodinger cat states. The 
nonclassicality 5 increases with the separation between the classical points defining the cat state. This growth 
saturates, if the separation distance is so large that the quantum interference patterns are well isolated from both main 
peaks of the distributions. Moreover, for a non-zero momentum po ^ 0, the quantity S undergoes oscillations until the 
separation distance becomes so large that both packets are separated from the interference patterns. Asymptotically, 
if the separation is large enough, the indicator 5 does not depend on po and tends to a constant value, 6max ~ 0.636. 

In the case of Fock states |n), the quantity 5 equals zero for the coherent vacuum state |0) and grows monotonically 
with the quantum number n. If a quantum state is displaced by the Glauber operator D{P), the shape of the Wigner 
function and the nonclassicality parameter do not change. Although the squeezing operator S{r]) changes the shape 
of the Wigner function, our results obtained for the squeezed Fock states show that the nonclassicality S does not 
depend on squeezing. 

The results presented in this work were obtained for pure states of infinite dimensional Hilbert space with use of 
the standard harmonic oscillator coherent states. It is worth to emphasize that our approach is also suited to analyze 
mixed quantum states. Furthermore, one may study the similar problem for quantum states of a finite dimensional 
Hilbert space, which was originally tackled in |37| . In such a case one defines the Husimi function with the help of 



the SU{2), spin coherent states, while the Wigner functions may be obtained by expanding the density matrix in the 
complete basis of the rotation operator s l4]l 14^. [4a| . The Wigner function for finite dimensional systems may also be 
defined in alternative ways - see ^S, HM HE lifl HE HE HEl and references therein. Studying the volume of the negative 
part of the Wigner function, defined according to any of these approaches, one may get an interesting information 
concerning the nonclassical properties of the state analyzed. For instance some recent attempts |36l ISll l52l | try to 
link the negativity of the Wigner function with the entanglement of analyzed quantum states defined on a composed 
Hilbert space, or with the violation of the Bell inequalities. 
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